TWO-PROTON CORRELATIONS 
FROM PB+PB CENTRAL COLLISIONS 

F. Wang (for the NA49 Collaboration) 

Nuclear Science Division, Lawrence Berkeley National Laboratory 

One Cyclotron Road, Berkeley, CA 94720, USA 

FQWang@lbl.gov 

Abstract The two-proton correlation function at midrapidity from Pb+Pb cen- 
tral collisions at 158 AGeV has been measured by the NA49 experiment. 
The preliminary results are compared to model predictions from proton 
source distributions of static thermal Gaussian sources and the trans- 
port models of rqmd and VENUS. We obtain an effective proton source 
size er c ff = 4.0 ± 0.15(stat.)^Q°g(syst.) fm. The rqmd model underpre- 
dicts the correlation function (a c s = 4.41 fm), while the VENUS model 
overpredicts the correlation function (cr e g = 3.55 fm). 
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Nuclear matter under extreme conditions of high energy density has 
been extensively studied through high energy heavy ion collisions. The 
baryon density plays an important role in the dynamical evolution of 
these collisions. To measure the baryon spatial density, one needs in- 
formation on the space-time extent of the baryon source. The space- 
time extent of the proton source at freeze-out can be inferred from two- 
proton correlation functions. The peak in the correlation function at 
Qmv = yZ—QfiQ^/^ ~ 20 MeV/c is inversely related to the effective size of 
the proton source [1, 2]. In the above, is the difference of the proton 
4-momenta, and qi nv is the momentum magnitude of one proton in the 
rest frame of the pair. 

We report the first, preliminary results on the two-proton correlation 
function in the midrapidity region from Pb+Pb central collisions at 158 
GeV per nucleon. The measurement was done by the NA49 experi- 
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ment [3] at the SPS, using the 5% most central events, corresponding to 
collisions with impact parameter b < 3.3 fm. 

Two independent analyses were performed on the data: dE/dx anal- 
ysis which used particles in the rapidity range 2.9 < y < 3.4 (assuming 
proton mass) with at least 70% probability to be a proton obtained 
from their ionization energies deposited in the time projection chambers 
(TPCs), and TOF analysis which used identified protons by combin- 
ing the time of flight information and the dE/dx. Both analyses used 
protons up to a transverse momentum px = 2 GeV/c. 

The proton samples are contaminated by weak decay protons (A + S° 
and S + ) which are incorrectly reconstructed as primary vertex tracks. 
From the measured single particle distributions [4, 5] and model calcula- 
tions of RQMD and venus, we estimate the contamination to be 15±5 5 %. 
This results in about 301^q% of the proton pairs having at least one 
proton from weak decays. These protons do not have correlation with 
protons from the primary interactions. In the dE/dx analysis, further 
contamination is present from kaons on the lower tail of their dE/dx 
distribution merging into the region where particles have at least 70% 
probability to a proton. This resulted in 25% K + p pairs and fewer than 
2% K + K + pairs in the proton pair sample. The coulomb hole of K + p 
pairs does not affect the resulting two-proton correlation function in the 
interested low qi nv region below 50 MeV/c. 

The two-proton correlation function is obtained as the ratio of the 
ginv distribution of true proton pairs to that of mixed-event pairs with 
protons from different events. The number of mixed-event pairs was 
large enough so that the statistical error on the correlation function is 
dominated by the statistical uncertainty in the number of true pairs. 
To eliminate the effect of close pair reconstruction inefficiency, a cut of 
2 cm was applied on the pair distance at the middle plane of the TPC 
for both true and mixed-event pairs. 

The correlation functions obtained from the two analyses can be di- 
rectly compared because of the nearly symmetric acceptances about 
midrapidity y c . m . = 2.9. The correlation functions (with the K + p con- 
tamination corrected in the dE/dx analysis) are consistent. In the results 
reported below, the true pairs and the respective mixed-event pairs from 
the two analyses were combined. The combined sample had about 10 5 
pairs with qi nv < 120 MeV/c, 75% of which were from the TOF analy- 
sis. The gi nv distributions of the true pairs and the mixed-event pairs 
are shown in the top panel of Fig. 1.1. The number of mixed-event pairs 
was normalized to that of true pairs in the range qw > 500 MeV/c. The 
resulting correlation function is shown in the bottom panel of Fig. 1.1. 
The prominent peak at q- mv ps 20 MeV/c is evident. There is a statisti- 
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cally significant structure in the correlation function at (fc nv rj 70 MeV/c. 
Many systematic effects have been studied; none has been identified that 
can account for the structure. There have been suggestions that a sharp 
edge in the two-proton density distribution of the source can produce 
such an effect [6]. 
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Figure 1.1 Top: q inv dis- 
tributions of true proton 
pairs (points) and mixed- 
event proton pairs (his- 
togram). Bottom: the 
measured two-proton cor- 
relation function. The con- 
tamination from weak de- 
cay protons and the finite 
momentum resolution were 
not corrected. The errors 
shown are statistical only. 
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We correct the measured correlation function for the contamination 
from weak decay protons and for the finite momentum resolution. The 
effect of the momentum resolution is only significant in the first two data 
points of the measured correlation function. The corrected correlation 
function is plotted in Fig. 1.2 as filled points. 

In order to assess the proton freeze-out conditions, we compare the 
measured two- proton correlation function to theoretical calculations. 
Given the proton phase space density distribution, the two-proton corre- 
lation function can be calculated by the Koonin-Pratt Formalism [1, 7]. 
Two types of proton source were used: 

(I) Isotropic Gaussian sources of widths a XtV ,a z and at for the space 
and time coordinates of protons in the source rest frame, and ther- 
mal momentum distribution of temperature T. No correlation between 
space-time and momentum of the protons is present. Following com- 



4 



binations of parameters were used in the calculations: a X;V = a z = a, 
a t = and a, and T = 120 MeV (as derived in [8]), 300 MeV (mea- 
sured inverse slope of proton transverse mass spectrum [4]) and 70 MeV 
(inverse slope observed at low energy, as an extreme). 

(II) Protons generated for Pb+Pb central collisions (b < 3.3 fm) at 
158 AGeV by two microscopic transport models: the Relativistic Quan- 
tum Molecular Dynamics (rqmd) model (version 2.3) [9] and the VENUS 
model (version 4.12) [10]. Both models describe a variety of experi- 
mental data on single particle distributions reasonably well. Protons at 
freeze-out have correlations between space-time and momentum intrinsic 
to the dynamical evolution in the models. 

Only protons in the experimental acceptance are used to calculate 
the two-proton correlation functions, the results of which are shown in 
Fig. 1.2 for rqmd, venus, and the Gaussian source with a x>y = a z = 
a t = 3.8 fm and T = 120 MeV, respectively. While venus overpredicts 
the amplitude of the correlation function, rqmd slightly underpredicts 
the amplitude. 



Figure 1.2 The two- 
proton correlation function 
after corrections for the 
30% contamination due 
to weak decay protons 
and the finite momentum 
resolution (points), com- 
pared to calculations for a 
Gaussian source (dotted), 
and for freeze-out protons 
from rqmd v2.3 (solid) 
and VENUS v4.12 (dashed). 
The errors shown on the 
data points are statistical 
only. 



We use x 2 / n df, the normalized mean square of the point-to-point 
difference between the data and the calculation in the range gi nv < 
48 MeV/c (i.e., 8 data points or ndf = 8), to quantify how well the 
calculations agree with data. We characterize the effective size of the 
proton source from each model by a e g = tyaAx • &Ay • VAz/V^, where 
cr Ax, & Ay and a& z are the Gaussian widths fitted to the distributions 
in Ax, Ay and Az, distance between the protons of close pairs with 
(frnv < 48 MeV/c. The distance is evaluated in the pair rest frame at the 
time when the later particle freezes out [11]. Respectively for rqmd and 
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venus, the x 2 / n df values are 1.47 and 2.29; the fitted Gaussian widths 
are (a Ax , a Ay , a Az ) = (5.91, 6.00, 6.83) fm and (4.57, 4.57, 6.08) fm, where 
z is the longitudinal coordinate; the resulting effective sizes are a c s = 
4.41 fm and 3.55 fm. 

In Fig. 1.3, we study the x 2 / n df as a function of a c g. The x 2 / n df 
values from all three models follow roughly the same solid line, drawn 
through the points for the Gaussian sources with T = 120 MeV to 
guide the eye [12]. From the minimum x 2 /ndf point and the points 
where x 2 / n df has increased by 0.125 (note ndf = 8), we extract a c g = 
(4.0 ± 0.15) fm, where 0.15 fm is the statistical error [13]. By applying 
a correction to the measured correlation function using a proton pair 
contamination of 20% and 50%, we obtain a systematic error of toig fm 
on <7 e ff. We note that u e g = 4.0 fm corresponds to a uniform density 
hard sphere of radius \fh(T c Q = 8.9 fm. 



Figure 1.3 The xVndf 
values as function of the 
effective source size a e s for 
various model calculations 
with respect to the mea- 
sured correlation function. 
The model calculations 
are for Gaussian sources 
(circles, triangles, and 
stars) and for freeze-out 
protons generated by 
RQMD v2.3 (filled square) 
and VENUS v4.12 (open 
square) . 



We have also studied Gaussian sources with extreme shapes: oblate 
&x,y = 4cr 2 and prolate a x>y = a z /4 (both with at = 0). The correspond- 
ing x 2 /ndf versus a e g curves are shown in Fig. 1.3 as dashed and dotted 
lines, respectively. They do not fall along the solid line, implying that 
the two-proton correlation function has certain sensitivity to the shape 
of the proton source. In principle, multi-dimensional two-proton corre- 
lation functions could reveal the shape of the proton source. However, 
the lack of statistics has prevented us from drawing conclusions. 

We note that there is no simple relation between <7 e ff and the pro- 
ton source in RQMD and venus model. From single proton distribu- 
tions at freeze-out, we obtain the following Gaussian widths in the 
source rest frame: (a x , a y , a z , at) = (7.6,7.7,6.4,7.0) fm for rqmd and 
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(3.6,3.6,4.3,1.9) fm for VENUS. The two-proton correlation function, 
therefore, appears to measure a smaller region of the source, i.e., a/\i < 
= x,y,z), due to space-time-momentum correlation. The effect 
is more dramatic in rqmd than in venus, which is consistent with the 
expectation that more secondary particle interactions in RQMD result in 
a stronger correlation between space-time and momentum of freeze-out 
protons. The fact that the x 2 /n&i values versus a c g for rqmd and VENUS 
lie on the curves obtained for the isotropic Gaussian sources, in which 
no space-time- momentum correlation is present, suggests that the effect 
of the space-time- momentum correlation is small in a c s . 

Finally, we comment on our two-proton correlation function in the 
context of other measurements. The pion source size measured by in- 
terferometry increases with the pion multiplicity [14], which increases 
steadily with bombarding energy in similar colliding systems [15]. Due to 
the large pion-nucleon cross-section, one would expect that protons and 
pions freeze-out under similar conditions, therefore, the proton source 
size would increase with bombarding energy as well. However, our mea- 
surement, in conjunction with preliminary results obtained at GSI [16] 
and AGS [17] energies, shows that the peak height is rather insensi- 
tive to the bombarding energy. This implies that the effective sizes of 
the freeze-out proton sources are similar in heavy ion collisions over a 
wide energy range. More detailed studies are needed to understand the 
possible acceptance and instrumental effects in these measurements. 

In summary, the NA49 experiment has measured the two-proton cor- 
relation function at midrapidity from Pb+Pb central collisions at 158 
GeV per nucleon. From comparisons between the data and the cal- 
culations, we extract an effective proton source size of <r e gf = 4.0 ± 
0.15(stat.)^Q l j ) g(syst.) fm. The rqmd model underpredicts the ampli- 
tude of the correlation function (<7 e ff = 4.41 fm), while the VENUS model 
overpredicts the amplitude (a e g = 3.55 fm). Due to the space-time- 
momentum correlation, the two-proton correlation function is sensitive 
only to a limited region of the proton source. Our measurement together 
with the measurements at lower energies suggest a very weak dependence 
of the two-proton correlation function on bombarding energy. 
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